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DOMAIN WALL ROUGHENING IN DISORDERED MEDIA: 
FROM LOCAL SPIN DYNAMICS TO A 
CONTINUUM DESCRIPTION OF THE INTERFACE 
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QWe study the kinetic roughening of a driven domain wall between spin-up and spin-down domains 
for a model with non-conserved order parameter and quenched disorder. To understand the scaling 
£Nj ■ behavior of this interface we construct an equation of motion and study it theoretically. 
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I. INTRODUCTION 



Roughening phenomena of growing surfaces and moving interfaces are of great interest for more than one decade 
since they appear in a variety of technological and scientific problems. The most interesting topic in the field of surface 
growth is the molecular beam epitaxy (MBE) due to its technical importance [Q. Examples for the roughening of 
moving interfaces are the immiscible-fluid displacement in oil recovery and the dynamics of magnetic domain walls, 
see for example [^J. Typical for the interface problem is the occurrence of the so-called depinning transition ||. 
Depending on a driving force the interface gets trapped by random impurities, i.e. it is in the pinning phase. For 
values above the critical force the interface moves steadily with non-zero velocity, it is in the depinning phase. In the 
examples listed these phenomena occur on very different length scales. For a description of all of these phenomena 
the concepts of self-affine fractals and dynamic scaling are of crucial importance [Q. 

Generally it is supposed that the surfaces as well as the interfaces can be described by an equation of motion 
t-H ! (EOM) which belongs to one of the famous universality classes, namely the Edwards- Wilkinson (EW) |5j or the 
■ Kardar-Parisi-Zhang (KPZ) j(3| universality class. In the first case the EOM for the surface (interface) profile function 
| h(x,t) reads 

O dh 2 

(N ■ — = i/V 2 h + r? + F (1) 



' with a noise term r\ and the driving force F. In the second case the so-called KPZ nonlinearity -|(Vh) 2 with A 
proportional to the surface (interface) velocity occurs additionally in Eq. (|l]). In the case of growing surfaces the 
■ noise arises mainly due to statistical fluctuations of the deposition rates and thus it can be assumed that the noise is 
time dependent (annealed) . In contrast the noise in the other examples listed arises from the moving of the interface 
I ■ through a random background, represented by a time independent (quenched) noise term. 

In the present article we study the motion and the morphology of a domain wall in a ferromagnetic medium 
containing quenched random fields. A continuum description of the interface is given and possible equations of 
^ ' motion arc discussed. 

> 

X 

b 

To study the behavior of a domain wall in a ferromagnetic medium (low temperature phase) we start from a 
Ginzburg-Landau type Hamiltonian 

H = J d d v (-|#r, tf + ^(r, tf + 3 - (V0(r, t)f - (H + B(r)) #r, *)) (2) 

where <f> denotes a scalar order parameter, H denotes a homogeneous driving field and B(r) a quenched random field. 
This continuum description of a ferromagnetic material on a mesoscopic scale is well established. It is assumed that 
a = b = uq without loss of generality and that the random fields, drawn with equal probability, have zero mean 
and are uncorrelated in space. Furthermore we assume J > corresponding to ferromagnetic spin-spin coupling and 



II. MODEL 
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wo > for stability. The dynamics of the system for non-conserved order parameter is defined through a Langevin 
equation 



dfat) on 

7 3t d<t>{r,t) [ ' 

with a relaxation time proportional to 7 (model A in the classification of Hohenberg and Halperin [Q), for further 
details see ^,^|. Temperature can be neglected since from renormalization group studies it is believed to be irrel- 
evant flo|| . Another argument for the irrelevance of temperature has been given by Grinstein and Ma |llf| for the 
random field Ising model (RFIM) in d dimensions using scaling arguments for the interface roughness w (defined in 
the next section). They found that it scales for temperature T — as w 2 ~ L 2 ( 4_:D )/ 3 with domain size L in the 
interface dimension T> = d — 1. In the pure Ising model it scales as L 2 ~ v for T > 0. Since 2(4 — 2?)/3 is for all T> 
greater than 2 — T> the field-induced perturbations dominate the thermal perturbations and thus temperature can be 
neglected. This argument can be applied to the present model since for a 00, b — > 00 and a/b = 1 one recovers the 
RFIM ||. 



III. DYNAMIC SCALING OF MOVING INTERFACES 



We assume a geometry where a horizontal domain wall separates one ferromagnetic domain with positive magneti- 
zation below the wall and one with negative magnetization above the wall. The position y — h (x, t) of the interface is 
then defined as that point r = (x, y) at which 0(x, y, t) as function of y for fixed x changes sign. Numerical studies || 
starting from a discretized version of Eq. (|^) showed that h (x, t) is a single valued function for not too large driving 
and random fields. 

Starting from an initially flat interface the interface develops a rough structure due to the random fields which 
increase or decrease the driving force locally. This roughness can be measured by the root-mean-square fluctuations 
of the averaged interface position Q 

u;(i,<) = ((h(x,i)-(h(x,0)) 2 ) 1/2 (4) 

which is usually called the roughness. Here L denotes the system size perpendicular to the moving direction of the 
interface and the angular brackets denote an average over all lattice sites at positions x as well as over different 
realizations of the disorder. Another quantity of interest is the height correlation function 

CM) = <[h(x + r,i)-h(x,i)] 2 >, (5) 

which is related to the roughness w(L, t) according to 

w(L,t) 2 = ^J2C(r,t). (6) 

r 

Eq. (^) is exact for periodic boundary conditions and for open boundary conditions for L — > 00. We will concentrate 
in the following on the scaling behavior of the height correlation function as a characteristic function describing the 
behavior of the interface. 

Due to the interaction terms in the underlying model the single interface elements do not grow independently of 
their neighbors and it is assumed Q that there exist a time dependent correlation length parallel to the interface 
direction. For r <C the correlations of the interface do not grow further and their saturation value is assumed to 
scale like 

C(r«£(t),i)cxr 2Q (7) 

with the roughness exponent a. For r larger than the time dependent correlation length the interface fluctuations 
are uncorrelated in space but they are still growing with time, their time evolution can be described by 

<7(r»£(i),*) oct 2 ^. (8) 

These scaling relations (0) and (||) are considered as limiting cases of a general dynamic scaling relation 

CM) = (£(i)) 2Q 3 (^) . (9) 
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In order that Eqs. (f7j) and (||) are recovered the scaling function g(y) must have the properties g(y) f» const for y 3> 1 
and g(y) oc y 2a for i/<Cl. Furthermore the correlation length has to grow with time according to 

oc i 1/z (10) 

with the dynamic exponent z — a/ (3. At time r ~ L z £(t) is of the order of the system size and it cannot grow further 
thus a crossover to a finite size behavior occurs. 

From these scaling assumptions follows that the interface scales as 

showing that the interface can be conceived as a statistical self-amne fractal on length scales x <C Here h(x, i) = 
h(x, t) — (h(x, t)) denotes the deviation of the interface from its time dependent average and f(p) denotes a scaling 
function with the limiting properties \ f(p)\ ~ const for /;> 1 and /(p) ~ p Q for p — > 0. 



IV. CONTINUUM DESCRIPTION OF THE INTERFACE 



We start the discussion of Eq. (||) by considering first the case B(r) = 0. The corresponding Langevin equation (j^) 
reads 

= _ Uo ^ (r; t ) (0( rj i) 2 - l) + J A 0(r, i) + JT . (12) 

Due to a driving field H > the domain wall moves in the positive y-direction and an approximate ansatz for the 
solution of Eq. (g) is 

0(r,t)=V (» + «*) (13) 
where translational invariance perpendicular to the y-direction is assumed. The velocity v of the domain wall is given 

by mi 

v s ^ff^o (14) 

where f2o = \/2J/uq is the intrinsic width normal to the domain wall in the static problem = (iJ = 0). 

Taking local curvature elements of the domain wall into account an ansatz 

K*,t)=*(*^p2) (15) 



n 

where h (x, t) denotes the time dependent wall position and 

n = fioVl + (Vh) 2 (16) 
the local width of the wall in y-direction has been assumed p^| . Within this ansatz a KPZ-like EOM 

|J = z,V 2 h + (Vh) 2 + n (x, h(x, *)) + F (17) 

with v = v(J), effective magnetic fields r\ and F and the noise correlator 

(77(x , h )?7(xo + x, h + h)) = D5 v (x)5(h) (18) 

can be obtained for a medium with random fields, see Ref. (14) for a similar approach. For an analysis of the EOM 
it is reasonable to substitute the time derivative <9h / dt on the right side of Eq. (|l^) by the averaged interface velocity 
v so that Eq. ( ]l7| ) can be replaced by 

^ = i ,V 2 h + ^(Vh) 2 + ?7 (x ! h(x,£))+ J F (19) 

where A/2 denotes the product of A and v. 
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V. DISCUSSION OF THE EOM 



In the quenched disorder case considered here the exponents characterizing the scaling behavior are not known 
exactly. We will show that with special assumptions at the depinning transition, F ~ Fc, and for the fast moving 
interfaces, F ^> Fc, the values of the exponents can be obtained. 

At the depinning transition the average velocity of the domain wall is zero. Therefore Eq. (|l^) is reduced to the 
quenched form of the EW equation 

— = ^V 2 h + ?/ (x,h(x,i)) + F. (20) 

The average over the noise at the pinned location of the interface, (r)(x, h(x, t))), compensates the driving force F at 
the depinning transition so that only the noise fluctuations <5f?(x, h(x, t)) = 7](ic, h(x, t)) — (r)(x., h(x, t))) are relevant 
in the EOM which then reads 

^ = i/V 2 h+fy(x,h(x,f)) . (21) 

For the correlations of the noise fluctuations we again assume 

(<fy(x 0) h )<5r/(x + x, h + h)> = D<P(x)5(h) . (22) 

Describing h as a generalized self-affine function with time dependent correlation length, as discussed above, a 
rescaling transformation 

h(x,t) = a~ a h(ax,a z t) (23) 

with a scaling factor a > 1 can be used to obtain the roughness exponent a and the dynamic exponent z. Under such 
a scaling transformation the noise, according to Eq. (^2j) scales as 

5f)(x, h) = a (p+Q)/2 ^(ax, a Q h) . (24) 



Inserting now Eq. (E3|) and Eq. (|24j) into Eq. (21) the EOM is recovered with rescaled coefficients 



= a*~ 2 v, D = a 2z - v - 3a D. (25) 

With z — 2 and a = (2z — f)/3 the coefficients v and D remain invariant and this invariance of the coefficients 
determines the values of the exponents. Therefore, the small time exponent [3 has the value f3 — a/z — (2z — T>)/3z. 
In simulations of the spin model |@ Eq. (@) we found near the depinning transition a w 0.9, z w 2 in P = 1 and 
a w 0.68, z w 1.89 in T> = 2 in good agreement with our theoretical values obtained from the scaling analysis. 

Our exponent relation for the roughness exponent a is in agreement with the prediction of Grinstein and Ma H ] 
who found for the RFIM a = (4 — T>)/3 (see above) . On the other hand, Leschorn |M found for an automaton model 



of Eq. (20) in the interface dimension V = 1 the values a « 1.25 and (3 ~ 0.88 and a » 0.75 and (3 « 0.48 in T> = 2. 
Renormalization group (RG) studies of Eq. ( pp| ) by Nattermann ct al. [p| in P = 4 — e lead to a = e/3 in agreement 
with our calculation and 2 = 2 — 2e/9 which disagrees with our ^-independent value for the dynamic exponent. 
Leschhorn notes that the dynamic exponent z obtained numerically from a measurement of a and (3 through z — a/ (3 
is very close to the RG values for different V. With dimensional analysis methods Parisi |l9fl a lso found j3 — (4 — T>)/6 
but this relation was questioned by him due to numerical results of an EOM similar to Eq. (|20D presented in a previous 
paper ^0|. He also gave arguments in favor of the scaling relation (3 — (4 — D)/4 in Ref. (19). In contrast, Kessler et 
al. HJ] as well as Dong et al. [B2] found numerically a » 1 for V — 1 and Ji et al. ^3] a « 0.67 for T> = 2 close to our 
values. However, discrepancies between the results of Dong et al. and of Ji et al. and our results still exist. Measuring 
the velocity dependence of F — Fc Dong et al. found v oc (F — Fc) 6 with 9 fa 0.25 while we found earlier || 9 « 1. 
Ji et al. found no self-affine scaling in V = 1. 

Since our results, both numerically and theoretically, for the roughness exponent a agree with other numerical work 
but disagree for the small time exponent (3 we also integrated numerically the EOM Eq. (GO) in T> = 1 in order to 
understand this discrepancy. We integrate a system of size L = 4096 up to a time t « 3.3 lO^vithin an Euler-scheme 
with dt = 0.02 and as the lattice spacing in h-direction we chose dh — 0.001. For the noise rj uniformly distributed 
in [—0.2,0.2] we observed a depinning transition close to F = 0.012. The results presented here are obtained for 
this driving force. Fig. (la) shows the spatial correlations of the height correlation function C(r,t) as function of the 
distance r for different times on logarithmic scales. An r-independent saturation is observed for large r in agreement 
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with Eq. (g) and a linear dependence on ln(r) for small r in agreement with Eq. (^). For the corresponding exponent 
we found a = 0.98 in good agreement with our theoretical result. However, the crucial point is that ln(C(r, t)) is still 
i-dependent for small r, i.e. contrary to the scaling prediction Eq. (R) the correlations are still growing with time, 
C(l,t) oc t 2K with k = 0.2 (see Fig. (lb)). This was also found previously in simulations of the corresponding spin 
model H Eq. (^). Note that such a behavior can be found in many different roughening models p3j. One possibility 
to take this anomalous behavior into account is a modification of the scaling relation Eq. (||) in the following way || 
(see Fig. (lc) for a scaling plot): 

C(r,t)=Z(t) 2a/X 9^) , (26) 

with the limiting forms 

C(r»at),t)^at) 2a/X (27) 

and 

C(r < f (i), t) OC r 2«^(f)2a(l/A-l) (28) 

where we have set A = 1 — k/(3 for convenience. Note that the exponent a again follows from the slope of ln(C) versus 
ln(r). 
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FIG. 1. Simulation results of Eq. (|2C|) for system parameters as mentioned in the main text, (a) The height correlation 
function C(r,t) for different times t, t increasing from bottom to top. (b) The time dependence of C(L/2,t) and C(l,£), the 
lines represent C(L/2,t) oc t 136 and C(l,t) oc t 0,4 . (c) Scaling plot of the height correlation function according to Eq. (p6[). 
(d) The time dependence of the correlation length £(t), the line represents a fit according to £(t) oc t l ' z with z = 2.18. 

If we define the dynamical exponent z as it is usual done in dynamic scaling theory by £(t) oc t x l z we have to set 
z = a//3X = a/(f3 — k). Therefore the exponent is no more the ratio of a and z. For the small time exponent (3 
we find from an analysis of the time dependent saturation value of the height correlation function at r = L/2 which 
according to Eq. (|27|) can be written as 
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C(L/2,t)«.t > * (29) 

from Fig. (lb) as well as from the width (not shown here) (3 = 0.68. From z = a/($ — k) we get z — 2.03 in good 
agreement with our theoretical result. 

The important result from this analysis is that only in the saturation limit in which (7(1, i) is time independent, 
i.e. on very large time scales and therefore for very large system sizes, the exponent (3 which describes the small time 
behavior of the height correlation function and of the width is given by the ratio a/z and the usual scaling behavior 
Eq. ([)]) is recovered. If C{\,t) is still growing proportional to t 2K the small time behavior (r <C £(i)) is described 
by the exponent (3 = n + a/z. Thus the discrepancy with the work of Parisi p0| mentioned above presumably has 
its origin in this modified scaling relation since with the time dependence of the width only the exponent (3 can be 
measured. 

To further support our point of view we additionally analyze the correlation length £(t) and thus the exponent z 
directly by the following approximative method. In a logarithmic plot of C (r, t) the quantity £ (t) can be obtained by 
the intersection of a horizontal line which is defined by the time dependent saturation value of the height correlation 
function and the line which one obtains by fitting C(r,t) for small r. With this method we obtain for the dynamic 
exponent z = 2.18 (see Fig. (Id)) which is of the order of the value we obtained from the relation z = a/((3 — k). 

The situation described so far is only valid at or very close to the depinning transition. For a moving interface with 
larger velocity v the KPZ nonlinearity A/2(Vh) 2 is important for the scaling behavior. For a fast moving interface, 
F ^> Fc, it is plausible that the quenched disorder acts as an effective annealed disorder so that in this limit the 
EOM is given by the well-known KPZ equation 

^ = vV 2 h + ±(Vhf + V (x,t) (30) 

where the constant driving force F was eliminated by going in the comoving frame of reference and the noise has a 
Gaussian distribution with zero mean and is uncorrelated in space and time 

{r)(x, t)r](x',t')) = D5{x - x')5(t - t') . (31) 

In this case the exponents are known exactly for T> = 1, namely a K pz = 1/2 and z K pz =3/2 thus /3 K pz = 1/3. In the 
annealed EW case these exponents are also known exactly and are given by a EW = 1/2 and /3 EW = 1/4 (z EW = 2). 
By simulation of Eq. ( |30| ) one observes a crossover at a time (™~* KPZ from an EW regime with (3 = /3 EW to (3 = /3 K pz 
depending on the coefficients of the KPZ equation |p5| . Comparing the formulas for the width of the interface in the 
KPZ case, 



< PZ = c 2 



|A|i 



2/3 



(32) 



to the corresponding width in the EW case for an infinite system, 



^ w = ^--^^)(2^- (33) 
2tt v z EW — 1 

one obtains for the crossover time between an EW behavior for small times to a KPZ behavior for large times with 
the numerical estimate ci « 0.4 

t w^KP Z w 252is 5 D- 2 \\\- 4 . (34) 

This relation limits the possibility of observing KPZ exponents since (3 can be observed only for times t where 
the correlation length is much smaller than the system size. For the EW model the correlation length can be 
calculated exactly pBfl 

£(*) = {2vt) l ' z ™ . (35) 

Therefore, systems of length L where periodic boundary conditions are assumed (the maximum value of the correlation 
length is L/2) saturate at times t = L z ™ /8i>. Comparing this result with Eq. ( |34| ) the system size Lq which one 
needs to observe the KPZ regime has to be much larger than 

L ^ 45 W' (36) 
We think that such a crossover from EW to KPZ behavior also occurs in V — 2 where for the EW equation logarithmic 
dependencies of the characteristic functions are known exactly and algebraic behavior in the KPZ case is observed. 
In our earlier simulations ]8|,[l6| for V = 1 and V = 2 we observed EW behavior for all possible system sizes which we 
believe are due to these finite size effects. 
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VI. CONCLUSION 



To summarize, we argued that the motion of a domain wall in a ferromagnetic medium described by a Ginzburg- 
Landau type energy can be mapped onto the quenched form of the EW equation at the depinning transition and onto 
the annealed KPZ equation for driving fields much larger than the critical field. In the latter case the roughness and 
the dynamic exponent are known exactly for V = 1 but as discussed above they can only be observed for very large 
systems while for smaller systems the annealed EW behavior is observed. At the depinning transition we obtained 
the exponents from a scaling analysis and obtained values in agreement with our numerical results presented earlier 
in V = 1 and V = 2. 
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